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Search for New Physics at = 7 TeV in Hadronic Final States with 
Missing Transverse Energy and Heavy Flavor 
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A search for supersymmetric particles in events with large missing transverse momentum, heavy flavor jet candi- 
dates and no leptons {e,fi) in y/s = 7 TeV proton-proton collisions is presented. In a data sample corresponding 
to an integrated luminosity of 0.83 fb~^ recorded by the ATLAS experiment at the Large Hadron Collider, no 
significant excess is observed with respect to the prediction for Standard Model processes. Model-independent 
production cross section upper limits are provided in the context of simplified models as well as conventional 
limits. 



1. Introduction 

! Supersymmetry (SUSY) [3 is one of the most compelling theories to describe physics beyond the Standard 
Model (SM). In the framework of a generic i?-parity conserving minimal supersymmetric extension of the SM, 
the MSSM [2j, SUSY particles are produced in pairs and the lightest supersymmetric particle (LSP) is stable. 
In a large variety of models, the LSP is the lightest neutralino, xl, which is weakly interacting and is a possible 

[ candidate for dark matter. The coloured superpartners of quarks and gluons, the squarks (g) and gluinos (g), 
are expected to be copiously produced via the strong interaction at the Large Hadron Collider (LHC). The 
partners of the right-handed and left-handed quarks, and q^, can mix to form two mass eigenstates. These 
mixing effects are proportional to the corresponding fermion masses and therefore become important for the 
third generation. In particular, large mixing can yield sbottom (b-^) and stop (ti) mass eigenstates that are 
significantly lighter than other squarks. Consequently, bi and ti could be produced with large cross sections at 
the LHC, either via direct pair production or, if kinematically allowed, through gg production with subsequent 
g — > bib or g — > tit decays. Depending on the SUSY particle mass spectrum, the cascade decays of gluino- 
mediatcd and pair-produced sbottoms or stops result in complex final states consisting of missing transverse 
momentum (its magnitude is referred to as E!^^^^ in the following) and several jets, among which 6-quark jets 
(&-jets) are expected. 

] SUSY is searched for in final states involving i?™"**^, energetic jets, of which at least one must be identified as 
a &-jet and no isolated leptons (e or fi). The search is based on pp collision data at a centre-of-mass energy of 

] 7 TeV recorded by the ATLAS experiment Q at the LHC in 2011. The total data set included in the analysis 
corresponds to an integrated luminosity of 0.83 fb~^. 

I Two phenomenological MSSM scenarios are considered where the first and second generation squark masses 
are set above 2 TeV. In the first scenario, the bi is the lightest squark, to- > to- > to ~o , and the branching ratio 

for g — )■ bib decays is 100%. Sbottoms are produced via gluino-mediated processes or via direct pair production 
I and they are assumed to decay exclusively via bi where m~o is fixed at 60 GeV. The interpretation of 

the results is presented as a function of the gluino and light sbottom masses. The second MSSM-like scenario 
is defined in the context of the general simplified models 01 : all squarks including bi are heavy, gluino-pair 
production is the only kinematically allowed process and gluinos decay (off-shell) into bbxi final states. Here 
the results are interpreted in a (jUg, m^^o) plane. These results are generalized to any new physics process where 
gluino-like particles decay into bb and a weakly interacting massive particle. 



2. Monte Carlo Simulated Samples 

Simulated event samples are used to determine the detector acceptance, the reconstruction efficiencies and the 
expected event yields for signal and background processes. Samples of SUSY signal processes were simulated for 
various models using the HERWIG++ [91 v2.4.2 Monte Carlo program. The particle mass spectra and decay modes 
were determined using the SUSYHIT [6] vl.3 program, and parametrized in the {iti-, m- ) and {rrig, m- ) planes. 

The SUSY sample yields are then normalized to the expectations from next-to-leading order (NLO) calculations 
obtained using the PROSPINO ^] v2.1 program. For these calculations the CTEQ6.6M ^] parametrisation of the 
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Physics process 


(7- BR [nb] 


W ^ Iv (+jets) 
Z/^* U (+jets) 
Z ^ vu (+jets) 
ti 

Single top 


31.4±1.6 [9-nJ 
3.20±0.16 [9-nJ 
5.82±0.29 [9-nj 
0. 1651° °} J [12-i^ 
0.085±0.003 [15,^6] 



Table I: The most important background processes and their predicted cross sections, multiplied by the relevant branching 
ratios (BR) before any event selection. A generator level cut mu > 40 GeV was applied to the Z/7*(— >■ t^l') process. 
Contributions from higher order QCD corrections are included for W and Z boson production (NNLO corrections) and 
for tt production (NLO+NNLL corrections). 



parton density functions (PDFs) is used and the renormalization and factorisation scales are set to the average 
mass of the sparticles produced in the hard interaction. 

For the background, the foUovifing Standard Model processes are considered: 



tt and single top production: events were generated using the generator MCONLO [17|, ll8| v3.41. For 
the evaluation of systematic uncertainties, additional ti samples were generated using the POWHEG [19.] . 



ALPGEN \M and ACERMC 21 



programs. 



W{-^ ^i/)+jets, Z/^*{-^ £+£~)+jets (where £ = e, /i, r) and 2(— > vv) +jets production: events with 
light and heavy (6,c) flavor jets were generated using the ALPGEN v2.13 program. A generator level cut 
mu > 40 GeV was applied to the Z/7*(— >• process. 

Di-boson {WW, WZ and ZZ) production: events were generated using ALPGEN, however, compared to 
the other backgrounds their contribution was found to be negligible, after the application of the selection 
criteria. 

For the QCD background, no reliable prediction can be obtained from a leading order Monte Carlo simulation 
and a data-driven method is used to determine the contribution to the selected event samples, as discussed in 
Section [31 

All signal and background samples were generated at y/s = 7 TeV u sing the ATLAS MClO parameter 
tune [22], processed with the GEANT4 23] simulation of the ATLAS detector |24l | . then reconstructed and passed 
through the same analysis chain as the data. For all generators, except for PYTHIA, the HERWIG + JIMMY [sl. [25j 
modelling of the parton shower and underlying event was used (v6.510 and v4.31, respectively). 

For the comparison to data, all non-QCD background cross sections are normalized to the results of higher 
order QCD calculations. A summary of the relevant cross sections is given in Table ID For the next-to-next- 
to- leading order (NNLO) W and Z/j* production cross sections, an uncertainty of ±5% is assumed. For 
the ti production cross section, the corresponding uncertainty on the NLO+NNLL (next-to-next-to-leading 
logarithms) cross section is estimated to be 

All Monte Carlo samples are generated with both in-time and out-of-time pile-up from multiple proton-proton 
interactions. The simulated events are reweighted such that the distribution of interactions per crossing in the 
Monte Carlo matches the one observed in data. 



3. Data and Baseline Event Selection 

After the application of beam, detector and data-quality requirements, the data set used for this analysis 
correspond to a total integrated luminosity of 0.83 fb"^. 

Events are selected at the trigger level by requiring one jet with high and large missing transverse 
momentum. The selection is fully efficient for events containing at least one jet with pT > 130 GeV and 
Sf""" > 130 GeV H^. 

In the data sample selected, jet candidates are reconstructed using the anti-kt jet clustering algorithm [27l - l29j 
with a distance parameter of i?=0.4. The inputs to this algorithm are three dimensional calorimeter energy 
clusters. The reconstructed jet energies are corrected for inhomogeneities and for the non-compensating nature 
of the calorimeter by using p^- and ?7-dependent calibration factors determined from Monte Carlo simulation 
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and validated using extensive test-beam measurements and studies of pp collision data (Ref. [30| and references 
therein). Only jets with px > 20 GeV and within jryj < 2.8 are retained for this analysis. Candidates for 6-jets 
are identified among jets with px > 50 GeV using an algorithm that reconstructs a vertex from all tracks which 
are displaced from the primary vertex and associated with the jet. The parameters of the algorithm are chosen 
such that a tagging efficiency of 50% (1%) is achieved for 6-jets (light flavor or gluon jets) in ti events in Monte 
Carlo simulation [3lj. 

In order to ap ply the lepton veto, electron candidates are required to satisfy the 'medium' selection criteria, as 
detailed in Ref. [32l |. Muon candidates are identified either as a match between an extrapolated inner detector 
track and one or more segments in the muon spectrometer, or by associating an inner detector track to a muon 
spectrometer track. The combined track parameters are derived from a statistical combination of the two sets 
of track parameters. Electrons (muons) are required to have px > 20 GeV (10 GeV) and jryj < 2.47(2.4). 

The calculation of iJ™"^*' is based on the modulus of the vectorial sum of the px of the reconstructed jets 
(with px > 20 GeV and over the full calorimeter coverage |?7| < 4.9), leptons (including non-isolated muons) 
and the calorimeter clusters not belonging to reconstructed objects. 

After object identification, overlaps are resolved. Any jet within a distance Ai? — 0.2 of a medium electron 
candidate is discarded and any remaining lepton within Ai? = 0.4 of a jet is discarded. 

Events are selected if the primary vertex is associated with five or more tracks. They are also required 
to pass basic quality criteria to discriminate against detector noise and non-collision backgrounds. Due to a 
front-end electronics failure in one of the electromagnetic calorimeter modules, a region of the calorimeter of 
size A77 X A0 = 1.6 x 0.4 was only partially read out. Events with any jet with px > 50 GeV in this region are 
rejected. The acceptance loss caused by this selection cut is about 10%. 

Selected events are required to have at least one jet with px > 130 GeV, at least two additional jets with 
Px > 50 GeV and E^^^^ > 130 GeV. At least one jet is required to be 6-tagged. Events containing electron or 
muon candidates are rejected. The effective mass, mes, is defined as the scalar sum of E!^^^^ and the transverse 
momenta of the three leading jets. Events are required to have E^^^'^/m.eg > 0.25. In addition, the smallest 
azimuthal separation between the E^^^^ direction and the three leading jets, A0„Mn, is required to be larger 
than 0.4. The last requirement effectively reduces the amount of QCD background where E^^^^ results from 
mis-reconstructed jets or from neutrinos emitted along the direction of the jet axis by heavy flavor decays. 

4. Signal Region Optimization 

Phenomenological models have many advantages as signal samples when attempting to design analyses sensi- 
tive to as broad a slice of phase space as is both possible and practical. They describe well-motivated, SM-like 
production and decay processes, and the kinematics are determined by a small number of parameters (masses). 
Both SUSY scenarios considered in this analysis result in 4 6-jet + E^^^^ final state signatures. The simplified 
model kinematics are determined to first order by a single parameter: the mass difference between the gluino and 
the neutralino, AM(g — xl). This simplicity motivates the choice of this model for signal region optimization 
studies, with cross-checks performed to ensure the results are relevant in the more complicated case. 

The optimization procedure was designed to ensure broad sensitivity given the baseline selection constraints. 
Several likely kinematic variables were evaluated for signal/background separation power: E^^^^, mes, jet 
multiplicity, b-tagged jet multiplicity, etc. From this, the best variable set was chosen and an n-dimensional 
cut grid was created. For each set of cuts and simulated signal sample, a systematic-corrected significance was 
computed (systematics estimated from Monte Carlo) and used to assemble a set of optimal cuts for maximal 
sensitivity to the entire (rng,m^o) plane. 

This list of "optimal" signal regions included nearly unique sets of cuts for each mass point. The number of 
selections were gradually reduced while ensuring broad-based sensitivity was retained. Four signal regions were 
chosen to represent good compromises among mass plane coverage, sensitivity, and practical concerns such as 
background control regions. They are characterized by the minimum number of &-jets required in the final state 
and by the threshold of the further selection on mgff: 3JA (>1 fe-jet, moff > 500 GeV), 3JB (>1 6-jet, moff 
> 700 GeV), 3JC (>2 6-jet, meff > 500 GeV) and 3JD (>2 6-jet, meff > 700 GeV). 

5. Standard Model Background Estimation 

The expected amount of tt, W/Z +ieta and single top events is estimated using the Monte Carlo simulation. 
Events from ti production represent the largest background component in all four signal regions. The Monte 
Carlo prediction is validated by a data-driven estimate which relies on control regions with the same kinematic 
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selection on jets and missing transverse momentum and an electron or a muon with > 20 GeV in the final 
state, meff > 600 GeV and 40 GeV< mx <100 GeV (where mx is the transverse mass computed from the lepton 
4-vector and the E^^^^) and at least one or two 6-jets. The normalization determined in these control regions 
(corrected for non-tt contamination) is then transferred to the kinematically similar signal regions. Figure [T] 
shows the mgff distribution in the 1-electron and 1-muon (1 6-tag) control regions. The agreement between the 
data and the Monte Carlo prediction is good and, as a consequence, the prediction of the data-driven estimation 
agrees with that of the Monte Carlo. 
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Figure 1: Effective mass distribution for the 1-electron (left) and 1-muon (right) tt control regions. The lower plot shows 
the bin-by-bin ratio of the data to the MC. The yellow band shows the full systematic uncertainty on the SM expectation. 

The Monte Carlo estimation of the W/Z background, which is done using ALPGEN as baseline generator, 
includes a dedicated simulation of W+ heavy flavor quarks. Double counting arising both at generator and at 
Matrix Element/Parton Shower matching level between samples generated with light and heavy flavor quarks 
is resolved using a AR matching between partons and jets and a dedicated overall normalization scale factor is 
applied to the W+ heavy flavor samples. For each signal region, the normalization of the inclusive W/Z Monte 
Carlo prediction is validated with a combined fit of tt and W/Z background components to the distribution of 
the number of 6-tagged jets in a 0-lepton control region defined by reverting the mcff cut. The fit confirms the 
Monte Carlo prediction. 

Since its contribution to the total background is small, the estimation for the single top background is based 
entirely on the Monte Carlo prediction. 

The remai ning QCD background in the signal regions is estimated with a data driven procedure. The 
technique [H, Hjj used is to smear the momentum of jets in clean data events with low E^^^^ to generate "pseu- 
doevents" with possibly large E^^^^ values. The method was validated by comparing data and pseudoevents 
distributions in QCD enriched control regions that are kinematically similar to the signal regions, obtained by 
reverting the cut on A(j)min- Figure [2] shows the mes distribution for events with A0„iin < 0.4 and 1 (left) or 2 
(right) 6-tagged jets. The emulated QCD distributions agree with the data. 



6. Systematic Uncertainties 

The background from top and W/Z production is obtained using the Monte Carlo simulation. The total 
uncertainty on this prediction is estimated to be between ±30% and ±35% depending on the final selection. 
It is dominated by the uncertainty on the jet energy scale, on the theoretical prediction of the background 
processes and on the determination of the 6-tagging efficiency. The uncertainty on the jet energy scale (derived 
using 2010 collision data "30^) varies as a function of the jet pt and pseudorapidity and it is about 2% at 
Pt = 50 GeV in the central detector region. Additional systematic uncertainties arise from the dependence of 
the jet response on the number of expected interactions per bunch crossing and on the jet fiavor. The total jet 
energy scale uncertainty at 50 GeV in the central detector region is about 5%. This translates into a 20-25% 
uncertainty on the absolute prediction of the background from SM processes. Uncertainties on the theoretical 
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Figure 2: Effective mass distribution for the QCD control region (A0mi„ < 0.4) for events with 1 (left) or 2 (right) b- 
tagged jets. The lower plot shows the bin-by-bin ratio of the data to the SM expectation. The QCD prediction validated 
in these control samples is derived using the smearing method described in the text, other SM background contributions 
are estimated from MC. The yellow band shows the full systematic uncertainty on the SM expectation 



cross sections of the background processes (see Section [5]) , on the modelling of initial and final-state soft gluon 
radiation and the limited knowledge of the PDFs of the proton lead to uncertainties of ±25% and ±30% on 
the absolute predictions of the tt and the W/Z+jet backgrounds, respectively. An additional uncertainty of 
50%(100%) is assigned to the associated production of W{Z) and heavy flavor jets. The uncertainty on the 
determination of the tagging efficiency for 6-jets, c-jets and light flavor jets introduces further uncertainties on 
the predicted background contributions at the level of ±10% (±22%) for tt and ±15% (±30%) for VF/Z+jets 
in the 1 (2) 6-tag signal regions. For the QCD background estimation, the uncertainty of 50% is dominated by 
the dependency of the smearing function on the flavor composition of the low E^^^^ sample used as smearing 
starting point. 

For the SUSY signal processes, various sources of uncertainties affect the theoretical NLO cross sections. 
Variations of the renormalization and factorisation scales by a factor of two result in uncertainties of ±16% for 
gg production and ±30% for b^b^ pair production, with little dependence on the sparticle masses and the SUSY 
model. 

The number of predicted signal events is also affected by the PDF uncertainties, estimated using the 
CTEQ6.6M PDF error eigenvector sets at the 90% C.L. limit, rescaled to la. The relative uncertainties on the 
gg {bib^) cross sections were estimated to be in the range from ±11% to ±25% (±7% to ±16%) for the gg {bibi) 
processes, depending on the gluino and sbottom masses. Uncertainties due to the modelling of initial and final 
state radiation on the signal are not included. 

The impact of detector-related uncertainties, such as the jet energy scale ( JES) and 6-tagging uncertainties, on 
the signal event yields depends on the masses of the produced sparticles. The total uncertainty varies between 
±20% (±35%) and ±10% (±10%) for the 1 (2) &-tag case as the gluino/sbottom masses increase from 200 GeV 
to 1 TeV, across the different scenarios, and it is dominated by the JES and the 5-tagging uncertainty for low 
and high mass sparticles, respectively. 

Finally, an additional ±4.5% uncertainty on the quoted total integrated luminosity was taken into account, 
based on the 2010 luminosity calibration which was transferred to the 2011 data by using the LAr forward 
calorimeter and the tile calorimeter current measurements. 



7. Results 

Good agreement between data and Monte Carlo prediction is observed in the distributions of the moff, the 
E^^^^ and the pt of the leading jet as shown in Figure |3] before the moff cut for the signal regions with one and 
two 6-tags. 

The observed and predicted event yields in the four signal regions are given in Table |TT] together with the total 
estimated uncertainty on the predictions. The value used for the top background estimate includes the tt and the 
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single top contributions as predicted by the Monte Carlo. The ti component is validated using a dedicated data 
driven procedure. The M^/Z+jets background is estimated using Monte Carlo, and the uncertainties correspond 
to those discussed in Section El QCD-multijet contributions are estimated with the jet smearing method. The 
SM predictions agree with the observed number of events in all four signal regions. 



Sig. Reg. 


Data (0.83 fb"^) 


Top 


W/Z 


QCD 


Total 


3JA (1 btag mcft >500 GeV) 


361 


991 +82 


121 ±61 


15 ±7 


356+^f 


3JB (1 btag mcff >700 GeV) 


63 


07+15 
■^'-12 


31 ±19 


1.9 ±0.9 




3JC (2 btag niefi >500 GeV) 


76 


55_22 


20± 12 


3.6 ±1.8 


7911 


3JD (2 btag iTicfi >700 GeV) 


12 


7.811.1 


5±4 


0.5 ±0.3 


13.01^1 



Table II: Summary observed and expected event yields in the four signal regions. The QCD prediction is based on the 
jet smearing method described in the text. Systematic uncertainties for the Standard Model predictions are given. 

Since no excess with respect to the SM predictions is observed in the data, the results are translated into 
95% C.L. upper limits on contributions from new physics. Limits are derived using the CLg [1^ method, while 
the power constrained limit (PCL) j37| method is used for comparison with previous ATLAS results. 

The results are interpreted in terms of 95% C.L. exclusion limits for the SUSY scenarios described in the 
introduction. In Figure 3] the observed and expected exclusion regions are shown in the (m~ , m- ) plane for 

the hypothesis that the lightest squark is produced via gluino-mediatcd or direct pair production and decays 
exclusively via — &x°. The NLO cross sections are calculated using PROSPINO. For each scenario, the signal 
region resulting in the best expected exclusion limit is used: the selection 3JD provides the best sensitivity in 
most cases. If AM{g — 6^) < 100 GeV, signal regions with 1 6-tag are preferred, due to the lower number of 
expected 6-jets above pt thresholds. The regions 3JA and 3JB provide the best sensitivity when ^ m~ 

and sbottom pair production dominates. All systematic uncertainties on the signal and background contribu- 
tions are taken into account in these limits and include the fully-correlated detector- type uncertainties (JES, 
^-tagging, trigger, pile-up effects, luminosity) as well as the theoretical uncertainties on the signal (renormal- 
ization/factorization scale and PDF). Gluino masses below 720 GeV are excluded for sbottom masses up to 
600 GeV. The exclusion is less stringent in the region with low AM{g — &i), where low E^^^^ is expected. This 
search extends the previous ATLAS exclusion limit in the same scenario by about 130 GeV [38:] (180 GeV if 
using the same limit setting procedure). 

Results are also interpreted in the context of simplified models. In this case, all the squarks are heavier than 
the gluino, which decays exclusively into three-body final states (bbxi) via an off-shell sbottom. Such a scenario 
can be considered complementary to the previous one. The exclusion limits obtained on the (m-,TO-o) plane 
are shown in Figure [5] for gluino masses above 200 GeV. For each combination of masses, the analysis providing 
the best expected limit is chosen. The selection 3JD leads to the best sensitivity for gluino masses above 400 
GeV and AM{g — Xi) > 100 GeV. At low AM{g — Xi), soft 6-jets spectra and low i?™"^'' are expected, giving 
higher sensitivity to the signal regions 3JA and 3JB are preferred. Low gluino mass scenarios present moderate 
m.cs and high 6-jet multiplicity, thus favouring signal region 3JC. Neutralino masses below 200-250 GeV are 
excluded for gluino masses in the range 200-660 GeV, if AM{g — x^) >100 GeV. 

The results can be generalized in terms of 95% C.L. upper cross section limits for gluino-like pair production 
processes with produced particles decaying into bbxi final states. The cross section upper limits versus the 
gluino and neutralino mass are also given in Figure [5] 

8. Conclusions 

An update on the search for supersymmetry in final states with missing transverse momentum, 6-jet candidates 
and no isolated leptons in proton-proton collisions at 7 TeV is presented. The results are based on data 
corresponding to an integrated luminosity of 0.83 fb~^ collected during 2011 by the ATLAS experiment at the 
LHC. Events with at least three energetic jets, large E^^^^ and at least one 6-tagged jet are selected in four 
signal regions based on the number of 6-tagged jets (>1 or >2 5-jets) and on the value of effective mass (>500 
or >700 GeV). The dominant Standard Model backgrounds are estimated from Monte Carlo simulation and 
are validated with data. 

No excess above the expectation from Standard Model processes is found. The results are used to exclude 
parameter regions in various iZ-parity conserving SUSY models. Under the assumption that the lightest squark 
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bi is produced via gluino- mediated processes or direct pair production and decays exclusively via 6i — bxi , and 
that m^o = 60 GeV, gluino masses below 720 GeV are excluded with 95% C.L. for sbottom masses up to 600 GeV 
using the CLg approach. This extends the previous (35 pb~^) ATLAS limits on gluino masses in the same 
scenario by about 130 GeV. Results are also interpreted in simplified models, where gluinos decay into heavy 
flavor final states (66xi) via an off-shell sbottom. In these scenarios exclusion limits in the (m~,m-o) plane are 

derived and Xi masses below 200-250 GeV are excluded for gluino masses below 660 GeV, if AM{g — x?) >100 
GeV. 95% C.L. upper cross section limits for gluino-like pair production processes with produced particles 
decaying into bbxi final states are also given. 
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Figure 3: Distributions of the effective mass, rucff (top), E^^^^ (middle) and the pT of the leading jet (bottom) for data 
and the expected SM processes in the 1 6-tag (left) and 2 6-tags (right) signal regions. The yellow band shows the 
full systematic uncertainty on the SM expectation. For illustration, the distributions of one reference SUSY signal are 
superimposed. The lower plot shows the bin-by-bin ratio of the data to the SM expectation. 
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Figure 4: Observed and expected 95% C.L. exclusion limits in the {mg,m^ ) plane. Also shown are the 68% and 99% 
C.L. expected exclusion curves. For each point in the plot, the signal region selection providing the best expected limit 
is chosen. The neutralino mass is set to 60 GeV. The result is compared to previous results from ATLAS and CDF 
searches which assume the same gluino-sbottom decays hypotheses. Exclusion limits from the CDF and DO experiments 
on direct sbottom pair production are also shown. 




Figure 5: 95% C.L. upper cross section limits in pb and observed and expected limit contours in the (mg,m-o) plane for 
gluino masses above 200 GeV. For each scenario, the signal region selection providing the best expected limit is chosen. 



